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POLYNOMIAL TERM STRUCTURE MODELS 


SI CHENG AND MICHAEL R. TEHRANCHI 
UNIVERSITY OF CAMBRIDGE 


Abstract. In this article, we explore a class of tractable interest rate models that have the 
property that the price of a zero-coupon bond can be expressed as a polynomial of a state 
diffusion process. Our results include a classification of all such time-homogeneous single¬ 
factor models in the spirit of Filipovic’s maximal degree theorem for exponential polynomial 
models, as well as an explicit characterisation of the set of feasible parameters in the case 
when the factor process is bounded. Extensions to time-inhomogeneous and multi-factor 
polynomial models are also considered. 


1. Introduction 

A time-homogeneous factor model of the risk-free interest rate term structure is one in 
which the time-t spot interest rate is of the form 

n = R{Zt) 

and the time-t price of a zero-coupon bond of maturity T is of the form 

Pt{T) = H{T-t,Zt) 

where i? : / —)■ M and H : M_|_ x / —)■ M are given non-random functions and Z = {Zt)t>Q 
is a given time-homogeneous Markov process with state space I modelling some underlying 
economic factor. For such a model to be sensible, the functions R and H and the process Z 
must be intimately connected. The goal if this paper is to study this connection when the 
function H{T — t, ■) is assumed to be a polynomial. 

The motivation for this study is classical. Recall that a sufficient condition for the bond 
market to have no arbitrage is that there exists an equivalent probability measure Q under 
which the discounted bond prices P(T), dehned by 

are local martingales for all T > 0. Furthermore, a sufficient condition to ensure that the 
discounted bond prices are local martingales is the union of the following two assumptions: 
hrstly, an analytic assumption that the function R is a classical solution of the partial 
differential equation 

(1) d^H = ^ aij - RH on R+x I 

l<i<d 
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for a given spacial domain J C and for given functions 6 : / —)■ and a : I ^ with 

boundary condition 

(2) H{0, z) = 1 for all z E I; 

and secondly, a probabilistic assumption that Z is a weak solution of the stochastic differ¬ 
ential equation 

(3) dZt = b{Zt)dt + a{Zt)dWt, 

taking values in I for all t > 0, where hh is a Brownian motion under the hxed pricing 
measure Q and a = aa~^. 

In principle, the above partial differential equation ([^ with boundary condition (|^ can be 
solved numerically whenever the functions 6, a and R are suitably well-behaved. However, 
resorting to a numerical method to solve the partial differential equation can be too slow 
to use practice for the purpose of calibrating the parameters of the model, and indeed, it 
obscures the relationship between the dynamics of the factor process and the resulting bond 
prices. Therefore, there has been considerable interest in developing tractable models, where 
the function H is of reasonably explicit form. 

Perhaps the two most famous tractable factor models are those of Vasicek mi and Cox, 
Ingersoll & Ross |1] . In these models the factor is scalar and identihed with the spot interest 
rate, so in the notation above, d = 1 and R{z) = z, while the functions b and a are affine 
and the function H is of the exponential affine form 

H{x,z) = 

It is easy to see that the partial differential equation Q reduces to a system of coupled 
Riccati ordinary differential equations for the functions Hq and hi and the boundary condition 
(|^ becomes ho(0) = /ii(0) = 0. Furthermore, it is well-known that in both cases the 
corresponding stochastic differential equation ([^ always has a unique local solution. While 
the local solution to the Vasicek stochastic differential equation is in fact the unique global 
solution, the situation with the Cox-Ingersoll-Ross stochastic differential is more delicate: 
for some values of the parameters, local solutions may explode in hnite time by hitting 
the boundary of the state space. Duffie & Kan studied exponential affine models where 
the factor process is of arbitrary dimension d > 1, hnding conditions under which the 
corresponding stochastic differential equation ([^ has a non-explosive solution. Subsequently, 
there has been a considerable body of research on the properties of these exponential affine 
models. A notable contribution to this literature is a general characterisation of exponential 
affine term structure models by Duffie, Filipovic & Schachermayer |6]. 

An exponential affine model can be considered a special case of the family of exponential 
quadratic models. An early example of a quadratic model was proposed by Longstaff na, 
and has since been developed and generalised by Jamshidian 1121. Leippold & Wu liai. and 
Chen, Filipovic & Poor [2] among others. 

One may wonder if there exist non-trivial exponential cubic (or higher degree) models. 
Filipovic answered this question in the negative, by showing that the maximal degree for 
exponential polynomial models is necessarily two. That is to say, the exponential quadratic 
models are indeed the most general class of exponential polynomial models. 

In this article, we consider a related class of models, in which the function H{x, ■) itself is a 
polynomial. For instance, in the case where the factor process is scalar-valued, the function 
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H is of the form 


H{x,z) = '^gk{x)z'^ 

k=0 

for n + 1 differentiable functions Qk : M+ —)■ M. One of the main results is a classihcation 
of all such models when H is assumed to satisfy a partial differential equation of the form 
of equation Q. It turns out that the functions b, a and R are necessarily polynomials 
of low degree and the functions {gk)k solve a system of coupled linear ordinary differential 
equations. In light of Filipovic’s maximal degree theorem for exponential polynomial models, 
it might come as a surprise the degree n is not constrained. We further hnd necessary 
and sufficient conditions on the model parameters such that the corresponding stochastic 
differential equation (|^ has a non-explosive solution. 

This work is inspired by the interest rate model of Siegel [TB]. He showed that for all 
integers d > 1 there exist an explicit affine functions R and explicit quadratic functions b, 
such that the partial differential equation ([^ has a solution H such that H{x, ■) is affine for 
all a; > 0. Note that in this case dz-zjH vanishes identically, and hence the function a = aa~^ 
need not be specihed to verify the partial differential equation. Furthermore, he showed that 
for a certain choice of a that the corresponding the stochastic differential equation (|^ has a 
non-explosive solution valued in the bounded state-space 


/ = {{zi,..., Zd) : > 0 for all i and 


We mention also the rational model of Brody & Hughston [T]. Working under the objective 
measure P, the state price density is modelled Vt = a{t) + (3{t)Mt where a and (3 are 
deterministic functions and M is a P-martingale. We show in section that the Brody- 
Hughston model is a special case of the time-inhomogeneous polynomial models considered 
here. 

Just as the Brody-Hughston model and the Siegel model described above, most of the 
polynomial models of this paper (but not all - see section 3.2) have the property that the 


spot interest rate is hounded. This stands in contrast to many familiar models, such as the 
Vasicek and Cox-Ingersoll-Ross models. Nevertheless, the range of the spot interest rate can 
be expressed easily in terms of the model parameters, and hence the range can be calibrated 
to any desired (hnite) width. 

Finally, a related work is that of Cuchiero, Keller-Ressel & Teichmann [B], who characterise 
a class of time-homogeneous Markov process Y with the property that the n-th (mixed) 
moments can be expressed as a polynomial of the initial point Yq of degree at most n. 
Indeed, consider the d = 1 case and let be the family of polynomials of degree at most n: 


Fn=\f-f(z) = Y,hz\ ft€ 


k=0 


They study the processes Y that have the property that for any degree n and any polynomial 
g & Fn, for all f > 0 there exists a polynomial h E Fn such that 


E[g{Y,)\Yo = y] = h{y). 
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In contrast, in this work we study processes Z that have the property that for a jixed degree 
n and a fixed function i?, for all f > 0 there exists a polynomial h = H{t, ■) G such that 

E[e-^oRiZs)ds^^^ _ 

In particular, their results do not imply ours, or vice versa. For further existence results 
for multi-dimensional polynomial preserving processes, consult the recent paper of Filipovic 
and Larsson uni. 

In the remainder of this article is arranged as follows. In section we present one of 
the main results, a classihcation of scalar time-homogeneous factor models which satisfy 
both the analytic assumption that bond prices satisfy a certain partial differential equation 
and the algebraic assumption that the bond prices can be expressed as a polynomial of the 
factor. In section we show that a polynomial model in which the spot interest rate is 
bounded from below necessarily has the property that the spot interest rate is also bounded 
from above; we also present an explicit example of a polynomial model in which the spot 
interest rate is unbounded from below. In addition, we provide a complete classihcation of 
models satisfying the probabilistic assumption that the corresponding stochastic differential 
equation has a non-explosive solution valued in a bounded interval. This section contains 
an easy-to-check formulation of Feller’s test of explosion for stochastic differential equations 
with polynomial coefficients, which might have independent interest. In section]^ we present 
a spectral representation of the bond prices in the context of a polynomial model. In section 

we consider a concrete example of this class of models. The Cox-Ingersoll-Ross model 
is recovered from this example in the limit of degree parameter n 'I' oo. The properties of 
this example are analysed and its parameters calibrated to US interest rate data in the case 
n = 2. Finally in section]^ we briehy discuss two extensions: a Hull-White-type extension 
where the coefficients are allowed to be time dependent, and the higher dimensional case 
d > 1. 


2. An algebraic result 

This section contains one of the main result of this paper, a classihcation of models that 
satisfy the analytic assumption that the pricing function H solves a particular partial dif¬ 
ferential equation, in addition to having the extra structural property that H{x,-) is a 
polynomial of hxed degree. To more clearly see the structure of the argument we consider 
only the time-homogeneous case with d = 1 in this section. The time-inhomogeneous and 
multi-dimensional cases are considered in section The following theorem is of a purely 
algebraic, rather than probabilistic, nature in the sense the factor process is not mentioned. 
Related probabilistic results are stated in section 

Theorem 2.1. Fix n > 1, and suppose that 

n 

(4) H{x,z) = '^gk{x)z^ 

k=0 

for 77, -|- 1 differentiable functions Qk : M+ —)■ M. Furthermore, assume there exists functions 
b, a and R such that 

d^H = bd,H + -RH on R+x I 


(5) 
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where I is a non-trivial interval. 

If the functions {gk)k o,re linearly independent, then the following holds true: 
Case n = 1. 

(A) R{z) = Rq + Riz and h{z) = bo + biz + 62 ^^ where = 62 - 

(B) {go,gi) solve the system of linear ordinary differential equations 

ilo = ~Rogo + bogi 

gi = ~Rigo + (^1 ~ Ro)gi- 


Case n > 2 . 

(A) R{z) = Rq + RiZ + R2Z^, b{z) = bo + biZ + b2z‘^ + b^^z^ and (t‘^{z) = Oq + aiz + a2z‘^ + 

+ 042:"^ where the coefficients are such that 

R 2 = f&3 = and Ri = nb2 + 

(B) {go,... ,gn) solves the system of linear ordinary differential equations 


(jk —gk -2 ( {k — 2)63 + ^- ^04 — R2 


+ 9k-l ( {k — 1)&2 + 


{k-l){k-2) 


as — Ri ] + gk \ kbi + 


k{k-l) 


02 — Ri 


I n i\7 k{k 1) \ (/c + 2)(/c + l) 

+ 9k+i ( {k + l)&o H -) + 9k+2 -- ao, 


where we interpret g _2 = 9-i = 9n+i = 9 n +2 = 0 . 


Before proceeding to the proof, note that Theorem |2.1| has an obvious converse: If H 
has the polynomial form of equation Q for a certain degree n and that the functions R, 
b and have the form given by condition (A) and the functions {gk)k solve the system of 
ordinary differential equations given by condition (B), then the function H solves the partial 
differential equation (|^. The proof is straightforward. 


Proof. Let 


Ak{z) = kb{z)z^-^ + ^ — ^\^{z)z'^-^ - R{z)z^ 


Equation (|^ holds if and only if the equation 

n n 

(6) E gk{x)z^ = '^gk{x)Ak{z) 


k=0 


k=0 


holds identically. 
For all m > 0, let 


= <^ / ; J ^ M : f{z) = fo,...,fm& 


fc =0 


be the set of polynomial functions of degree at most m. Since / is non-trivial, a polynomial 
/ G Fm uniquely specifies its coefficients fo,. .., fm- 

We first show that if equation (|^ holds then the functions Ak G for all k. To see this, use 
the assumed linear independence of the functions {gk)k to pick n-l-l points 0 < xo < ■ ■ ■ < Xn 
such that the (n -|- 1) x (n -I- 1) matrix { 9 i{xj))ij is invertible. By evaluating equation (|^ 
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at the points (xj)j and solve for the Ai[z), we see that Ai{z) is a linear combination of 
monomials z’^ of degree at most n. 

Case n = 1 . Note that 


R{z) = -Ao{z) 
b{z) = Ai{z) + zR{z). 


Since Aq and Ai are in Fi, i.e. are affine, then R is affine and b is quadratic. Letting 
b{z) = bo + biz + b2Z^ and R{z) = Rq + Riz the above system equation implies 62 = -Ri- 
Finally, the identity ([^ becomes 

go + giz = go{Ro + Riz) + gi{bo + (61 - zRo)z). 

Equating coefficients of z yields the necessity and sufficiency of the system of ordinary 
differential equations. 

Case n > 2 . Note that 


R{z) = -Ao{z) 
b{z) = Ai{z) + zR(z) 
cr^(z) = A2 (z) — 2 zb(z) + z^R{z). 


Since the functions A^ are polynomials, so are the functions i?, b, and a^. On the other hand 


An{z) = nb{z)z^-^ + 



”-2 - R{z)z^ 



and, since the term in brackets is a polynomial, we have 


(7) 



Similarly, since A^-i G and An -2 ^ Fn we have 



a\z)-R{z)z^ e F3 C F4 


a\z)-R{z)z^ e F4. 


Since 



nb{z)z + 



(n-2)(n-3) 

2 


a\z) - R{z)z^ 


- 2 


^(n - l)b{z)z + 


inclusions ([^, (|^ and ([^ together yield 


( 10 ) 
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Similarly, since 


zb{z) = ( nb{z)z + 




a{z) — R{z)z^ ) ~ “ l)b{z)z + 


— {n — l)a^ 

inclusions ([^, (|^ and ( 10 ) together yield 

( 11 ) beFs. 

Finally, inclusions ([^, (10) and 10 together yield 

R e Fa. 

Recall that An is of degree at most n. Now substituting R{z) = Y^'l^oRk^^, b{z) = 
= J2k=o^kZ^ into the dehnition of An, and setting the coefficient of 

to zero yields 

( 12 ) 


n+2 


, n(n - 1 ) 
nbs H---04 = Fa 


Similarly, equating to zero the coefficient of z"'^^ in the expansion of An yields 

n(n — 1 ) 

noa H--03 = Fi- 


Finally, equating to zero the coefficient of z^~^^ in the expansion of An-i yields 
(13) 


/ (n — l)(n — 2) 

(n — 1)63 H--- 04 — Fa 


Note that equations (12) and (13) together are equivalent to 

n(n — 1 ) 


Fa = -63 = 

2 2 


-04. 


Finally, substituting these expressions into equation (|^ and comparing the coefficients of 
the monomials z^ yields the system of ordinary differential equations for the functions for 
the functions (gk)k- D 


3. Some probabilistic results 

In this section we include some results related to the probabilistic assumption that a 
certain stochastic differential equation has a non-explosive solution. 

3.1. Bounded state space. In this subsection, we argue that there are good reasons to 
make the further assumption that the state space I of the factor process F in a polynomial 
model is bounded, at least in the one-dimensional case. 

Recall that we aim to model the price Pt{T) at time t of a zero-coupon bond of maturity 
T by the formula Pt{T) = H{T — t,Zt) where Zt is the economic factor at time t. Since 
the payout of the bond is its face value Pt(T) = 1, there are economic grounds to assume 
that the bond prices are bounded. Indeed, to avoid a buy-and-hold arbitrage, one must have 
Ft(T) > 0; furthermore, assuming the existence of a bank account continuously paying the 
spot interest rate r* and assuming that rt > —C for some constant F > 0, then no arbitrage 
would imply that PtiT) < . 
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Continuing with this argument, suppose that the function H{x, ■) is bounded on the state 
space / C for each a; > 0. Furthermore, suppose that H{x, ■) is a non-constant polynomial. 
If the dimension d = 1, then by an elementary fact of real analysis we can conclude that I 
is bounded. 

We summarise the mathematical content of the above argument in the following proposi¬ 
tion: 


Proposition 3.1. Let I C W be a non-trivial interval. Suppose that for eaeh z & I, there 
exists an I-valued weak solution Z of the stochastie differential equation 

dZt = b{Zt)dt + a{Zt)dWt, Zq = 

Furthermore, suppose that there exists a funetion R sueh that 

H{x,z) = = z] 

for all X >0. Assume that R has the property that there exists a eonstant C > 0 such that 

R{z) > —C for all z E I. 

If H{x, •) is a non-eonstant polynomial for some a; > 0 then the interval I is bounded. 

Remark 3.2. Notice that the assumption that the spot rate = R{Zt) is bounded from 
below implies that the bond prices are bounded from above. In the case of a one-dimensional 
polynomial model, boundedness of H{x, ■) further implies the boundedness of the state space 
I. Finally, the boundedness of I and the continuity of R together the that spot interest rate 
is also bounded from above. 

On the other hand, notice that the boundedness of H{x, •) does not imply the boundedness 
of I in the case of exponential polynomial models such as Cox-Ingersoll-Ross. 


Remark 3.3. Note that in higher dimensions, non-constant polynomials may be bounded on 
unbounded sets. For instance, consider the polynomial 


on the set 


h{zi,Z2) = Zi- Z2 


I = {{zi,Z2) : |2:i - 2 : 2 ! < !}• 


We now mention a pleasant and well-known probabilistic consequence of the boundedness 
assumption. We state it here in the d dimensional case. 


Proposition 3.4. Let I C be bounded. Suppose the funetion H E satisfies the 
partial differential equation 

d:,H = ^ hi ^ X] aij - RH onR+x I 

l<i<d 


with boundary condition 

H{0, z) = 1 for all z E I] 

and suppose that for each z E I, there exists an I-valued weak solution Z of the stochastic 
differential equation 

dZt = b{Zt)dt -\- a{Zt)dWt, Zq = 
where a = . If R is continuous then 

H{x,z) = = z] 



for all X > 0 and z E I. 

Proof. Fix z E I, and let Z solve the stochastic differential equation with Zq = z. Also 
£x a time horizon x > 0. As mentioned in the introduction, since H satisfies the partial 
differential equation, then by Ito’s formula we know that the process M = dehned 

by 

Mt = e- ^0 - t, Zt) 

is a local martingale. Assuming that the state space I is bounded and that Zt E I for all 
t > 0, then process M is bounded by a constant by the continuity of R and H. 

In particular, the bounded local martingale M is a true martingale by the dominated 
convergence theorem, and hence 

H{x,z) = Mo = E[M,] 

as desired. □ 


3.2. An unbounded example. The message of Proposition 3T is that the assumptions 
that the spot interest rate is bounded from below and that the bond prices are polynomials 
in a scalar factor together imply that the spot interest rate is bounded from above as well. 

We note that many popular interest rate models respect the assumption that the spot rate 
is bounded from below - an example is the Cox-Ingersoll-Ross model. However, one should 
mention that some models, such as Vasicek, place no such lower bounds. In particular, while 
having a lower bound may be a reasonable and desirable feature, but we must grant that 
this assumption is not universally enforced. 

Therefore, we briefly explore the one-dimensional polynomial model, where we drop the 
assumption that the function R is bounded below on I. We will see via the following example 
that we can no longer conclude that the state space / is bounded. 

Let the state space be / = (0, oo), and the degree be n = 2, the coefficient functions be 
given by 

bE) = -Ft 


a{z) = z'^, 

and the bond pricing function be 


R{z) = 


Note that 


with initial condition 


H{x, z) = 1 + xz + -{e^ — X — 1 )^^ 


dxH = bd,H + -ad^^H - RH 


so we are in the setting of Theorem 2.1 


H{0,z) = 1 

By Ito’s formula, we know that if there exists a 


dWtj Zq — z 


solution Z of the stochastic differential equation 
(14) dZt =-^Z^ dt + Zt 

then the process 

Mt = H{T - t, Zt)e-^0 




is a local martingale. 

However, in this case, it is straightforward to show that the unique solution of the sto¬ 


chastic differential equation (14) is given by the formula 

^gm-i/2 

'Zit = 


1 -f I /g 

Furthermore, the local martingale M is a true martingale thanks to the identity 

. 2 " 

E{e^o = E 


-I 






= H{x, z) 

which can be verihed by explicit calculation. Note that there is no contradiction with 


Proposition 3.1 because inizei Riz) = —oo in this case. 


Remark 3.5. We mention here an interesting (though a bit tangential) observation regarding 
the above example. It is easy to see that if Z satishes the stochastic differential equation 


(14) then the process Y = dehnes a local martingale. Indeed, it satishes the stochastic 


differential equation 

dYt = Y^ogiYt) dWu Ho = eh 

It is slightly less obvious that the process H is a strictly local martingale. See, for instance, 
the paper of Goodman m for further results. 


3.3. A form of Feller’s test. As argued in Section 3.1, there are economic reasons to 


consider polynomial models in which the factor process takes values in a bounded inter¬ 
val. Therefore, in this subsection we consider solutions to the scalar stochastic differential 
equation 

dZt = b{Zt)dt + a{Zt)dWt, 

which live in a bounded state space / = (^ min . Furthermore, in light of Theorem 2.1 

we assume that the coefficients b and cx^ are polynomials. 

To avoid trivial complications, we assume 

= 0 = and a{z) > 0 for 2 ;min < ^ < 

Note that the coefficient b is Lipschitz on the closed interval [zmin, 2 ^max], while the coefficient 
a is Lipschitz on any interval [zmin + 1/-^, ^max — 1/-^] for A > 1 large enough. Therefore, 
for every z G (zmin, ^ max ) the stochastic differential equation has a unique nested family of 
strong solutions (^t, 7 v)te[o,rjv] with Zq^^ = z, where 

Tn = inf{t > 0: Zt^ [z^nin + l/N, ^^ax - 1 /^^)}- 
The explosion time T is then dehned as 

T = supT/v. 

N 

We are interested in the case where the solution is non-explosive in the sense that T = oo 
almost surely. 

The classical necessary and sufficient conditions on the functions b and a is Feller’s test 
of explosion. See for in stan ce Chapter 5 of Karatzas and Shreve •s na book for an account. 
Motivated by Theorem 2.1 we adapt Feller’s test to the case where the functions b and 
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are polynomials. It is likely that the following result is well-known, but we were unable to 
locate a reference in the literature. 


Theorem 3.6. Let b and a be polynomials. Furthermore, assume 

a(^min) = 0 = a(2:max) and a{z) > 0 for z^m < z < z^s.^. 

Letting a = y/a, there exists a unique non-explosive strong solution Z of the stoehastic 
differential equation 

dZt = b{Zt)dt + a{Zt)dWt 

taking values in the interval [zmin) ^^max) if and only if 

26(^min) a (^min) ^ 0 ^ 26(^jnax) ® (^max)- 


Proof. Recall that Feller’s test is 

P(T = CX)) = 1 n(^min) = 00 = n(2;max), 
where Feller’s test function is dehned by 


v{x) = 


eh a(^) dy dz, for Z^in < x < Zma.^, 


f z = c J y=z ) 


where c = kiz yuir, -+- z^^x)- h is enough to consider the behaviour of v near x = z^i^, as the 
behaviour near x = z^a^x is analogous. 

By changing variables, we now study the cases where the integral 

p{z) 


^rr\\n ) — 


a{z)p'{z 


-dz 


is hnite or inhnite, where 

p[z) = / e-^y dy. 

^ ^min 

is the related to the scale function. Now, by assumption the functions a and b are polyno¬ 
mials, and hence near Zmin can be written as 

a{t Zaiin) = 

b{t + z„,in) = /3t^ + 

for constants a > 0 and /3 ^ 0 and for integers A, B > 0. Note that with this notation 

26(^min) a (^min) ~ /31{B=0} Cll|^= 0 }. 

Hence, we must show that v{zr r,\A = oo on 

{H > 0, R = 0, /3 > 0} U {H > 0, R > 0} U {H = 0, R = 0,2/3 > a} 
and that n( 0 ) < cx) on the complement 

{H > 0, R = 0, /3 < 0} U {H = 0, R > 0} U {H = 0, R = 0, 2/3 < a}. 


We have the calculation 



2b{.s) 

a{s) 


ds 


const -|- 0(t) if i3 > H -|- 1 

— ^ log t -I- const + 0{t) ii B = A 
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and hence 


P (t ^min) 


const (1 + 0{t)) 

(const + 0 (t)) 


and therefore 


pit + Zrnin) 


ait '2^min)P d“ 


aB>A+i 
iiB = A 

ifB < A-1 

iiB>A + l 
if B = A, 2(3 > a 
= l + Oit)) if B = A,2P <a 


^ H-^i + Oit)) 
oo 


oo if B < A — 1, (3 > f) 

^r^(l + 0(t)) if B <A-1,I3 <f). 


From this, we see that v{z m\A = oo precisely on 

{B>A + l,A>l}A{B = A,2l3>a}A{B = A>l,2l3<a} 

U {/I > S + 1, /3 > 0} U {A > S + 1 > 2, /? < 0} 

from which the conclnsion follows. 


□ 


3.4. A canonical parametrisation of scalar polynomial models. We are now in a 

position to snmmarise the above resnlts to characterise the range of admissible parameters 
for which there exists a scalar polynomial model in which the factor process takes valnes in 
a bonnded open interval I. Since we can replace the state variable Z with Z = p{Z) and 
the interest rate fnnction R with R = Ro where p is an affine transformation, there 
is no loss of generality in hxing the state space I to be any given interval. Therefore, to 
simplify some calcnlations, in this section we will set I = (—1,1) and will refer to this as the 
canonical state space in the seqnel. 

we are interested in necessary and snfficient conditions on the 


In light of Theorem 2.1 


parameters 6 o, • • •, ^ 3 , ^o,..., 04 where the stochastic differential eqnation 


dZt = hiZt)dt + aiZt)dWt 

has a non-explosive solntion valned in the open interval (— 1 , 1 ) where 

hiz) = 60 + biz + + hz^ 

cr'^iz) = Oo + aiz + a 2 Z^ -f 032;^ -1- 042;"^. 


In order to enforce the condition cr(—1) = 0 = cr(l) we rewrite cx^ as 

a'^iz) = (1 - z^)icQ + ciz + C2Z^). 


In order to describe the parameter space, we hrst let 
(15) C = {(cq. Cl, C 2 ) : Co + CiZ + C 2 Z^ > 0 for all — 1 < .^ < 1}. 

This set can be described more explicitly: 


Proposition 3.7. 

C = {cq > 0, —Co < C 2 < Co, |ci| < Co -|- C 2 } U {co > 0, C 2 > Co, |ci| < 2^C 0 C 2 }. 
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Proof. For a fixed triplet (cq, ci, C2) let c{z) = cq + ciz + C2z‘^. Suppose that c{z) > 0 for all 
— 1 < 2 ; < 1. Note that c(0) = Cq, so Cq > 0 necessarily. 

By continuity, we have c(—1) = Cq — Ci + C 2 > 0 and c(l) = Cq + Ci + C 2 >0. Hence 
C 2 > —Co and |ci| < Cq + C 2 necessarily. Note that 

f{z) > Co - |ci||2:| + C2^^ 

> Co — (cq + C2)|^| + C2^^ 

= (co - k|c2)(l - 1^1). 


From this we see that if C 2 < Co then c{z) > 0 whenever \z\ < 1, and hence the condition 
|ci| < Co + C 2 is also sufficient. 

So now suppose that C 2 > cq. Letting 


2:0 = 



we have 0 < 2:0 < 1 and 

c(2;o) = (ci + 2^/^)zo■ 

Hence the condition ci > — 2 ^coC 2 is necessary. By considering c(—zo) we see that the 
condition ci < 2 ^/coC 2 is also necessary. Finally, writing 




we see that the condition |ci| < 2 ^coC 2 is sufficient as well. 


□ 


We now consider the function b. Let 


( 16 ) 13 co,cuc 2 = {(^ 0 , bi, & 2 , bs) : 2b{-l) - a'(-l) > 0 > 26 ( 1 ) - a'(l)}. 

Again, this set can be described more explicitly: 

Proposition 3.8. 

^C0,C1,C2 = { 1^0 + 62 + Ci| < —(61 + 63 + Co + C2)} 

The proof is straightforward. 

We can now parametrise the general scalar, time-homogeneous polynomial model with 
factor process taking values in (—1,1). For each n > 2, we let 

( 17 ) Pn = {(-Ro, -Ri, i?2, 60, 61, 62, 63, Co, Cl, C2) : (co. Cl, C2) G C, (60, 61, 62, 63) G Hco,ci,C2) 

Ri = nb 2 - |n(n - l)ci, R 2 = {n - 1)63 = l{n - l){n - 2 ) 02 } 

To summarise, we have proven that if 

(Rq, Ri, R2, bo, 61, 62,63, Co, Cl, C2) G Vn 
then for every G (— 1 , 1 ) the stochastic differential equation 

dZf = (bo + biZt + b 2 Z^ + b^Z^)dt + (1 — z'^)(co + CiZ^ + C 2 Zf)dWt, Zo = z 
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has a unique strong non-explosive solution taking values in (—1,1). Furthermore, the solution 
has the property that 

n 

= z] = '^gk{x)z'^ 


where the functions g^,... ,gk solve the the system of linear ordinary differential equations 


A, {k- 2 )(k- 3 ) 

gk —gk-2 {{k — 2)63 -|-^-04 — R2 


+ 9 k-l ( {k — 1)&2 + 


{k-l){k- 2 ) 


0-3 ~ Ri ] + 9 k \ kbi + 


k{k-l) 


02 — R{ 


fn iNr k{k + l) \ {k + 2 ){k + l) 

+ 9 k+i ( (^ + 1)A H -^- Oi j + gk+2 - - - Oq, 


where we interpret g _2 = 9-i = 9n+i = 9 n +2 = 0, subject to the boundary condition 

fi'o(O) = l, 5 'i( 0 ) = ■ ■ ■ = gn{ 0 ) = 0 . 


4. A SPECTRAL REPRESENTATION 

We are in the setting of the scalar time-homogeneous polynomial model with the factor 
process taking values in a bounded open interval I. The coefficient functions {gk)k are the 
solution of the system of differential equations which can be equivalently described as follows. 
Let S = be the (n -|- 1) x (n -|- 1) matrix with entries 

^j+k,j jbk+1 "F “ 2 ®fc +2 Rk 

and where Rk = bk = Ok = 0 when k < 0 and Rk = bk+i = ak +2 = 0 when k > 2. For 
instance, when n > 4, the matrix has the form 

— Rq bo Oq 

—Ri bi — Rq 26q -|- cii Sclq 

— R 2 &2 — Ri 2b\ -|- 02 — Rq 36o T 3oi 600 

63 — i ?2 2&2 oq — Ri 3&1 -|- 3a2 — Rq 4&o + 6 ai 

263 -|- 04 — R 2 362 T 303 — Ri 4&4 -|- 602 — Rq 

Now letting 

G{x) = {go{x),gi{x),...,gn{x)y 

The system of differential equations becomes 

G = SG, 

and, in particular, the solution can be expressed as 

G{x) = e^^G(O), 

where the boundary condition for zero-coupon bond pricing is given by 

G(0) = ( 1 , 0 ,..., 0 )^. 

It turns out that the matrix S has a nice property: 

14 





Proposition 4.1. The eigenvalues Aq, ..., A„ of S are real and satisfy 

Ai < — inf R{z). 

z£l 

for all i. 


To prove Proposition 4.1, we first prove a result on the existence of an invariant measure 
which may have independent interest. 


Proposition 4.2. Suppose (co,Ci,C 2 ) is in the interior of the set C defined by equation (15) 
and that h is in the set Bco,ci,c 2 defined by equation (16) . Then there exists a positive, 


inteqrable function f satisfying the differential equation 


with boundary conditions 


bf = ^{afy 


lim a{z)f{z) = 0 = \ima{z)f{z), 

24,-1 2tl 


where 


b{z) = bo + biz + b 2 z‘^ + b^z^ 
a(z) = (1 — z‘^)[cq + CiZ + C 2 ^^). 


Remark 4.3. If the function / is normalised so that f{z)dz = 1, then / is the unique 
invariant density for the diffusion Z with drift b and volatility a. That is, if the initial 
condition Zq is distributed with density /, then Z^ has the same distribution for all f > 0 . 


Proof. Since (cq, Ci, C2) is in the interior of C, we may assume that either 

-Co < C 2 < Co and |ci| < co + C 2 


or 

C 2 > Co, |ci| < 21 /C 0 C 2 }. 

In either case, the function c(z) = Cq + Ciz + C 2 Z^ is strictly positive on the closed interval 
[—1,1]. In particular, the function a(z) = (1 — z^)c(z) has simple roots at z = —1 and z = 1 
with 

a'(-l) = 2c(-l) > 0 > a'(l) = -2c(l). 

Furthermore, since (&05 ^ 2 , ^ 3 ) £ '^(co,ci,c 2 )) we have 

6(-l) > 0 > 6(1). 


Now any positive solution to the differential equation is of form 


m = 


c 

a{z) 


pz 2i)(s) , 
gJO 0(3) 


for \z\ < 1 , where C > 0 is a constant. 


Following the proof of Theorem 3T we focus on the left-hand end point z 
show that such an / is integrable and t/(f — 1) —0 as t 0. Writing 


b{t - 1) = f + 0{t) 
a{t — 1) = at + Oiff) 


—1. We must 
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a routine calculation shows that 


fit - 1) = ^ 

a 


from which the conclusion follows. 
We can now prove Proposition |4.1 


□ 


Proof. We can and do assume that have the made the canonical choice of state space I = 
(—1,1) introduced in section 3.4 Since S varies continuously with the model parameters, 
there is no loss of generality to assume that the parameters are in the interior of the Vn 
dehned by equation ([T^. By Proposition 4.2 there exists an invariant density /. 

Consider the inner product on dehned by 


n n 

^—n A —n J—1 


A+j 


f{z)dz 


i=0 j=0 

= j PiAfiAfiAdz 


where" is the linear operator sending the vector 

P= (P0,---,Pn) 

to the n + 1 degree polynomial 


PA) = 


PkZ 


k=0 


The key observation is that 


Note that 


Sp = 2'^P' ~ ^P' 


(P, Sq) = - j [\ap'q' + Rpq]f dz 
= {Sp,q) 

by integration by parts, where we have used the boundary condition 

lim a{z)f{z) = 0 = \im a{z)f {z). 

24,-1 2tl 

In particular, we see that S is symmetric with respect to this inner product and hence all 
eigenvalues are real. The inequality 

(p, Sp) <- j Rp^f dz 


< — inf R(z)(p,p) 
- 1 < 2<1 


implies the claimed upper bound on the spectrum. 

16 


□ 







Remark 4.4. Of course, the eigenvalues of the matrix S are the zeros of the characteristic 
polynomial which has degree n + 1. It is well know that there exists an explicit formula, 
discovered by Ferrari in 1540, for the zeros of quartic polynomials, and hence the eigenvalues 
of S can be expressed in a closed formula in terms of the matrix entries when n < 3. In 
particular, in this case, the bond pricing function H can be written, at least in principle, 
explicitly in terms of the model parameters. 

When n > 4, there is little hope for explicit formulae for the function H in terms of 
the model parameters. However, note that the matrix is sparse, in the sense that there 
are at most hve non-zero matrix entries per row. In particular, the product of the matrix 
exponential and the vector (1,0 ,..., 0)^ can be computed efficiently, and hence the lack 
of explicit formulae is not necessarily a prohibitive disadvantage. 


The proof of Proposition 4.1 shows that when there exists an invariant density /, then 


= MS 


where M = {Mij)ij is the (n -|- 1) x (n -f 1) positive dehnite matrix with entries 


Mij = f{z)dz. 


Suppose that the n -|- 1 real eigenvalues of the matrix S are Aq, ..., A„. Then the matrix S 
has the spectral decomposition 


n 

S ^ ^ Aj Ui Vi 

i=0 

where Ui is the right-eigenvector and Vi the left-eigenvector associated to the eigenvalue Aj, 
scaled such that 


Vi%Lj 


1 li i = j 
0 if i ^ j. 


For convenience, we choose the normalisation 


uj Mui = 1, 


and note that the left- and right-eigenvectors are related by 

Vi = uJ M. 

Now given the ith right-eigenvector Ui we can form the polynomial Ui{z) = 
Note that Ui is an eigenfunction of a certain differential operator 

{b{z)d^ + \a{z)d:,^ - R{z)) Ui{z) = Ajhj(x), 

and that the ith left-eigenvector Vi is related to hj by the formula 


Vi,k = J z^Ui{z)f{z)dz 


In particular, the bond pricing function takes the form 


i=0 
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where the function Qi is the (at most) n degree polynomial 


Qi{z) = Ui{z) / Ui{s)f{s)ds 


That is to say, the bond price can be seen to be a linear combination of the bond prices 
arising from n + 1 models with constant interest rates r = —Xi, where the coefficients Qi of 
the combination depend on the factor process. Note that by setting x = 0 we have 


i=0 


SO it is tempting to think of the numbers {Qi{z))i as probabilities; however, in general 
Qi{z) < 0 for some i and z & I, so such an interpretation is not always valid. An example 
where Qi takes negative values can be found in section below. 

In the general case, where the parameters are such that no invariant density exists, the 
matrix S is not necessarily diagonalisable. In this case, the bond pricing formula must be 
modified to 

n 

(18) H{x, z) = Qi{x, 

i=0 


where now the weight functions Qi are polynomials in both x and ^ and can be computed 
from the Jordan decomposition. An example where the matrix S is not diagonalisable is 
discussed in Section 3^- though strictly speaking, the setting is slightly different there since 
the state space of the stochastic differential equation (14) is unbounded. 

One consequence of formula (18) is that the long maturity interest rate can calculated as 

-—\ogH{x^z) = — maxAj. 


lim 


X 


for all z ^ I, unless the coefficient Qi{x, z) of the maximum eigenvalue is identically zero. 


5. An example 

In this section we explore a concrete realisation of a polynomial model. The purpose of 
this account is as a proof of concept and is not intended as an endorsement of this particular 
model over others. 

In the general polynomial framework, the function i? : J —)• M, mapping the factor process 
to the spot interest rate, is a quadratic function. In the following example, we assume that 
R is affine. By an affine change of variables, we can and will take the spot rate itself as 
the factor process. Note that this choice of parametrisation differs from the canonical choice 
introduced in Section [3^ 

We consider here a model where the spot rate is the solution of the stochastic differential 
equation 

drt = a{/3 - rt)dt + fn) (7 - 

This model is inspired by the Cox-Ingersoll-Ross: the parameter (3 > 0 plays the role 
of a long time mean level, the parameter a > 0 controls the rate of mean-reversion, and 
the parameter 7 > 0 is related to the spot rate volatility, in the sense that the inhnitesimal 
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variance of drt is approximately '-y'^rtdt when r* is very small. Indeed, sending n ^ oo formally 
recovers the Cox-Ingersoll-Ross stochastic differential eqnation. 

However, nnlike the Cox-Ingersoll-Ross mo del, the interest rate stays within the bonnded 
interval I = (0, Applying Theorem 3.6, we see that the appropriate Feller conditions 

2 


are 


2a(3 > 7 ^ and a{'^n — 2/3) > 7 


,n 


n 


To fnrther simplify this discnssion, we set the degree parameter n = 2, so that the set of 
feasible parameters becomes 

0 < /5 < 7 < \/2a/3. 

The corresponding matrix S introdnced in Section takes the form 




0 a/3 0 

— 1 —a 2a/3 + 7 ^ 
0 —1 — 2(0 + 7 ) 


from which bond pricing fnnction can be evalnated efficiently by the formnla 

H{x,r) = (1, r, 0, 0)"''. 

Notice if the parameters are snch that 

0 < /3 < 7 < \/2a/3 

then the process is ergodic in the pricing measnre Q, and its invariant density is given by 
the stationary solntion of the corresponding Fokker-Planck partial differential eqnation: 

1 jr 2b(p) 


fir) oc 


a{r) 


gJ0 ct{p )2 


OC r 


C-i 


( 7 -r) 


-C-2g-e/(7-r)^ 


where 


C = ^ and 0 = 2a (^1 - ^ 


T 


7 


To conhrm that this model can prodnce sensible looking yield cnrves, the parameters 
a, /3 ,7 have been calibrated to US Treasnry constant matnrity rates freely available from the 
Federal Reserve Economic Data (FRED) website [ 8 ]. The 251 daily yield cnrve observations 
for matnrities 

Xi G {1/12, 3/12, 6/12,1, 2, 3,5, 7,10, 20, 30 years } 

and dates 

tj G { 19 Angnst 2015, ..., 18 Angnst 2016 } 
were £t to this model. The parameter valnes 

a = 0.248,/3 = 3.1%, 7 = 12.9% 
seem to be a (local) minimnm of the penalty fnnction 

^ ^ (hobs(^^Uj) (Xj, Tj )) 
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subject to Feller feasibility condition, where the model yield is Y{x,r) = —^\ogH{x,r) and 
the spot rate rj on day tj is linearly interpolated from the two shortest dated yields via the 
formula 

f'j = 2^obs(l/12, tj) — -yobs(3/12, tj). 

Figure shows the fitted yield curve compared to the observed yield curve on the first and 
last days of the sample period. It is no surprise that the fit is better for left-hand end, as 
there are more observations at short maturities. Although the fit is not extremely impressive, 
it is probably sufficiently good for the modest goal of illustrating the potential utility of the 
class of polynomial models. 

Recall that the bond pricing function has two representations 

H{x,r) = go{x) -h gi{x)r -h g 2 {x)r^ 

= Qo(r)e^°" + Qi(r)e^^" + Q2{r)e^^^ 

where the exponential polynomial functions {gk)k solve the system 


go = al3gi 

91 = -go - oigi + (2a/3 -f 7^)5f2 

92 = - 9 i - 2 (a + 7)5(2 


with initial condition 

5o(0) = l, 5 i( 0 )= 52 ( 0 )= 0 . 

and the quadratic functions (Qi)i satisfy the eigenequation 

^r (7 - rfQi(r) + a(/3 - r)Q'(r) - rQi(r) = AiQi(r) 

with the normalisation such that 

[ Qi{r)f{r)dr = f Qi{rff{r)dr, 

Jo Jo 

where (Aj)j are the eigenvalues of the matrix S. Note that Qi{r) < 0 for some values of r. 
The graphs of the coefficient functions {gk)k and {Qi)i are shown in Figures and 
See the recent thesis [3] of Cheng for other calibrated examples. 


6. Extensions 


In this section, we will extend theorem 2.1 in two different ways: namely allowing time 
dependency and allowing a multi-dimensional factor process. We focus on the algebraic 
result relating the analytic assumption that the bond pricing function satisfies a certain 
partial differential equation to the assumption that the bond prices are polynomial in the 
factor process. We omit a discussion of the probabilistic assumption of whether a certain 
stochastic differential equation has a non-explosive solution. 
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Figure 1. 



6.1. Hull— White-type extension. As usual, by incorporating time-dependent parameters, 
we can hope to have a better model calibration. We introduce time dependency both in the 
dynamics of the factor process (Zt)t>o and the coefficient functions (gk)k- As one may expect, 
we will establish a similar algebraic result in this case. 

To be clear, we now consider a factor process (Zt)t>o to be a non-explosive solution valued 
in a non-trivial interval / C M to the time-inhomogeneous stochastic differential equation 

dZt = b{t, Zt)dt + a{t, Zt)dWt. 
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Figure 2 . 



Figure 3 . 



The spot rate is modelled as r* = R(t, Zt) and the bond prices are dehned by 

n 

Pt(T) = 

k=0 
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where : A —)• M are smooth deterministic functions satisfying the boundary conditions: 


^ 7 o(T,T) = 1 

gk{T, T) = 0 for all 1 < k < n 
where A = {(f, T) : 0 < t < T}. 

By adding the t component, the analytic assumption is this case is that the functions gk 
satisfy the partial differential equation 

n n 

(19) '^dtgk(t,T)z’^ = '^gk(t,T)Ak(t,z) for all 0 < f < T, z e / 

k=0 k=0 

where 


Ak{t, z) = R{t, z)z^ — kb{t, z)z 


k-i k{k 1) 2 


z)z' 


k-2 


The algebraic result in this case becomes the following: 

Theorem 6.1. Suppose that for each f > 0, the functions {gkit, ■))k are linearly independent 
and satisfy equation (19). 

Case n = 1. 

(A) R{t, z) = Roit) + Ri{t)z and b{t, z) = bo{t) + bi{t)z + b 2 {t)z‘^ where Ri{t) = b 2 {t). 

(B) {go,gi) solves the system of linear ordinary differential equations 

—dtgo = —Rogo + &o5'i 
~dtgi = —Rigo + (bi — Ro)gi. 

Case n > 2. 

(A) R(t,z) = Ro(t) + Ri(t)z + R 2 (t)z^, b{t,z) = boit) + bi{t)z + b 2 (t)z‘^ + b^.iCjz'^ and 
cT^(z) = ao(t) + ai{t)z + a 2 (t)z‘^ + ay,{t)z'^ + a 4 ^{f)z* where the coefficients are such that 

R2{t) = Ri{t) = nb2{t) + 

(B) {go,... ,gn) solves the system of linear ordinary differential equations 

-dtgk =gk -2 ({k- 2)&3 + ———^04 - R 2 


in {k-l){k-2) k{k-l) 

+ fi'fc-i ( {k — 1)62 H-- ^3 ~ .Ri ) A gk\ kbi H---02 — Ro 


in iN 7 k[k 1) \ [k 2) [k 1) 

+ 9 k+i ( {k + 1)60 H -- fli ) + fi'fc+2 - - 0,0, 


where we interpret g _2 = g-i = fi'n+i = fi'n +2 = 0. 


The proof is essentially the same as that of Theorem 2.1 so is omitted. 
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6.2. Brody—Hughston rational model. In the paper [T] of Brody & Hughston, the fol¬ 
lowing rational model is discussed. Let M be a positive martingale under the objective 
measure P, and suppose Mq = 1. Set 


where a and (3 are positive, continuously differentiable, deterministic functions. The idea is 
that Id is a model for the state price density. Therefore, bond prices are given by the formula 

PtiT) = ^wF{yT\Pt) 

Vt 

^ a{T)+p{T)Mt 
a{t) + 

Note that the bond prices are a rational function of the random variable Mt-, giving the 
model its name. Furthermore, by setting a{t) +/d(f) = Po{t) for f > 0, this model can match 
the initial term structure of interest rates. 

On the other hand, notice that we can write the bond prices as 

Pt{T) = go{t, T) + gi{t, T)Zt 


where the coefficients are dehned by 




m 


m 


and where we let 




be the factor process. In particular, this is an affine factor model and hence should be 
described by Theorem 6T We now carry out the verihcation under the assumption that 

dMt = Mt)MtdBt 


where 5 is a P-Brownian motion and u is bounded. 

In this framework, we can dehne the spot rate as 

rt = -dTPt{T)\T=t 

a{t) -1- 
a{t) + 

= Ro(t) + Ri(t)Zt 


where 


Note that 


R^it) 


m 

m 


and Riit) 


${t)a{t) — a{t)/3{t) 

W) 


dVt = ^d(t) -l- dt + l3{t)dMt 

= -Vt{rtdt + XtdBt) 
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where A* is the market price of risk dehned by 

* a{t) + (3(t)Mt 

= iy{t,Mt){a{t)Zt - 1 ) 


Since the process {Xt)o<t<T is bounded, we can dehne the equivalent risk-neutral pricing 
measure Q by 


dP 


r 

eh 


Vs ds 


Vt 


_ p — ^fo Xtdt+fg \tdBt 


to recover the usual pricing formula 

Pt{T) = 

Finally, we consider the dynamics of the factor process Z = V~^. By Ito’s formula we 
have 


dZt — Zt[{rt -|- X‘^)dt -|- XtdBt] 

= ihit)Zt + b 2 {t)Zl)dt + ait, Zt)dWt 


where 


blit) = 

b2it) = 


Kt) 

/3{t) 

$it)ait) — d(f)/3(f) 


cr(f, z) = u [t, 


m 

1 — zait) 
z/3{t) 


iait)z — 1 )^ 


and where the process iWt)o<t<T dehned by 


Wt = Bt+ / Xsds 
Jo 

is a Q Brownian motion by Girsanov’s theorem. In particular, notice that the drift is qua¬ 
dratic in Z and b 2 it) = as predicted by Theorem 6.1, while the volatility is determined 
by the function v. 


6.3. Multi-dimensional factor process. In this subsection, we will extend the polynomial 
model framework by allowing both the factor process {Zt)t>o and the background Brownian 
motion iWt)t>o to be multi-dimensional. To be more specihc, let iWt)t>o be a H-dimensional 
Brownian motion. Let {Zt)t>Q be the factor process taking values in / C M'^, assuming to be 
the (non-explosive) solution of the stochastic differential equation 

dZt = biZt)dt + a{Zt)dWt 

for some deterministic functions 6 : / —)■ and a : I ^ We dehne the dihusion 

function a = aa~^, and note that the only role played by the parameter D is as the upper 
bound on the rank of the matrix a{z). 
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For k = (fci,..., fcrf) G and z = (zi,..., G W^, we define the monomial as follows: 

k_ ki ka 

z — Zi 

For each m > 0, we dehne the following set of indices 

Km = {fc G : ki + ... + kd < m}. 

With this notation, we consider models where the bond price is of the form 

P,(T) = 'MT - t)Zf 

k^Kfi 

where the functions satisfy the boundary conditions 

^,(0) = liffc = (0,...,0) 

QkiS^) = 0 otherwise 

The spot rate is modelled as r* = R{Zt) for some deterministic function i? : / —)■ M. 

The partial differential equation Q gives rise to the condition 

(20) E gk{x)z'^ = ^ gk{x)Ak{z) 

keKn keK„ 

holds for any x > 0 and z E I, where the functions Ak are dehned as 

^k{z) = bi{z)d;,.z’^ + ^ ^ aij{z)d^^^^z^ - R{z)z^. 

l<i<d i<i,j<d 

Finally we dehne the notation 

Fm=\f-.I^^,f{z)=J2fk^^ /fceMi 

I k&Km ) 

to be the family of polynomials in d variables of total degree less or equal to m. We will as¬ 
sume that the state space J C is non-trivial in the sense that / G Fm uniquely determines 
the coefficients {fk)k, or more precisely so that 

= 0 for all G / implies /fc = 0 for all k G Km- 

k&Km 

For the sake of brevity, we only consider the case n >2. 


Theorem 6.2. Suppose n > 2 and that the functions gk are linearly independent and satisfy 


equation (20). Then we must have R G F 2 , h G F 3 for all 1 < i < d and Uij G F 4 for all 
^ F i, j < d. Furthermore, the coefficients are constrained in such a way that Ak G for 
all k such that k such that n — l<ki + ... + kd<XL. 


The proof follows the same pattern as the proof of Theorem 2.1 with heavier notation to 
account for the extra dimensions. We include it here for completeness. 

Proof. First we show that the functions A^ G Fn are polynomials for all k G K^. Let 
N = ("'^^) be the cardinality of set K^. Since the functions gk are linearly independent, we 
can hnd N distinct points Xi,..., independent of z such that the matrix with i-th column 
formed by vector {gk{xi), k G Kn) is non-singular. Now £x any 2;, we can rewrite condition 
(20) as a set of N simultaneous linear equations with N unknowns Ak{z). Therefore the 
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solution exists and is unique and can be written as linear combinations of the monomials 
hence all of the Ak{z) are polynomials in d variables of total degree less or equal to n. 

In what follows, we will hnd it convenient to introduce the notation 

(a)j = ( 0 ,..., 0 , a, 0 ,..., 0 ) where a is the i-th component. 

(a, b)ij = ( 0 ,..., a,..., 6 ,..., 0 ) where a is the i-th. component and b is the j-th component. 

for certain Z^^-valued indices, where a, 6 G Z_|_. 

Since we must have Ak{z) G Fn for all k G Kn, we can conclude for any 1 < i,j < d 

^o(^) = ~R{^) ^ Fn 

A(^i)Xz) = bi{z) - ZiR{z) G Fn 

A{i,i)i^^{z) = bi{z)zj + bj{z)zi + aij{z) - ZiZjR{z) G 

Therefore we may conclude immediately that the functions R{z), bi, aij are polynomials. On 
the other hand 

^(n).(^) = nz^-\{z) + ^^ z^-^aii{z) - z^R{z) G Fn 

by cancelling the z^~‘^ factor, we may deduce that: 

(21) nzibi{z) + aii{z) - z^Riz) G F 2 

Similarly by considering v 4 („_ 2 ),, ., we get 

(22) (n - l)zibi{z) + — —— —aii{z) - ZiR{z) G F 3 

(23) (n - 2)zibi{z) + — —— —auiz) - z‘fR{z) G F 4 


(24) {n-l)ziZjbi{z) + zfbj{z) + 


2j. , (^-2)(n- 1) 


Subtracting (21) from (22) and subtracting 


Zjaii{z) + {n- l)ziaij{z) 
from (|2^ gives 


zfzjR{z) G F 3 


Zibi{z) + (n - l)aii{z) G F 3 
Zibi{z) + (n - 2)aii{z) G F 4 


Hence we get the required degree constraint on functions R, 6, a. For the remaining part of 
this theorem, observe that given the degree constraint, the functions Ak will automatically 
be in Fn as long as ki + ... + kd < n — 2. □ 
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